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The recent progress on unquenched QCD simulations is critically reviewed. After some discussions on problems 
and subtleties in unquenched simulations, hadron spectra obtained from both quenched and unquenched simula- 
tions are compared among various gauge and quark actions. It is found that the Edinburgh plot does not agree 
in the continuum limit between Wilson and KS quark actions even in quenched QCD. Dynamical quark effects 
on hadron spectra, in particular, on meson masses are then presented for Wilson-type quark actions. Finally 
dynamical quark effects on other quantities such as the topological susceptibility and the flavor-singlet meson 
mass are discussed. 



1. Introduction 

An ultimate goal of lattice QCD simulations is 
to numerically solve the dynamics of the strong 
interaction. Due to lack of enough computa- 
tional resources, however, a large part of sim- 
ulations so far have been performed within the 
quenched approximation. Recent results from 
high statistical simulations 0| identified errors as- 
sociated with the quenched approximation: pre- 
dicted hadron masses in the continuum limit, in 
particular those containing strange quarks, differ 
from corresponding experimental values as large 
as 10%. The effect of the quark determinant must 
remove these differences. 

In this talk I review recent progress in un- 
quenched QCD simulations. I will survey the 
quality of current unquenched simulations and 
to what extent dynamical quark effects can be 
seen, in the light of the precise results from the 
quenched simulation 

I first discuss several problems in unquenched 
simulations such as the auto-correlation and the 
finite size effect, which slowed down a lot un- 
quenched simulations in the past. Although these 
problems still exist, they are not a main ob- 
struction anymore. A brief consideration on the 
scale determination and a theoretical discussion 
on the problem of changing the quark mass in 
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unquenched QCD then follow. 

Using a large part of this report hadron spec- 
tra are analyzed for KS and Wilson quark ac- 
tions in parallel, to compare the quality of data. 
Unfortunately and unexpectedly I notice that 
they disagree in the continuum limit even for the 
quenched case, as we will see. 

In the main part the effect of dynamical quarks 
on hadron spectra is discussed in the case of the 
Wilson quark action. After the chiral extrapo- 
lation in the sea quark mass, the effect becomes 
most visible, in particular on the hyper-fine split- 
ting of mesons, which consequently removes a 
large part of the quenching error in strange meson 
spectra. 

In the last part, the string breaking, the topo- 
logical susceptibility and the flavor-singlet meson 
mass are discussed in connection with the dynam- 
ical quark effect. This year mostly positive results 
are reported on all these quantities. 

Some important subjects such as quark masses, 
weak-matrix elements and heavy quark physics in 
unquenched QCD will not be covered. Please see 
other reviews on these topics[|]-||. 

2. Recent unquenched QCD simulations 

Recent large-scale unquenched simulations are 
summarized in Table |^, where Nf is the num- 
ber of flavors of dynamical quarks. Throughout 
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Table 1 

Summary of recent unqucnchcd QCD simulations 



action 



a(fm) 



La (fm) mps/mv 



Nf = 2 
PWpl 

pwjnf 

PKSllI 

PKs|l|] 

PC(1.76)J 

RC(TP)f 

PC(NP)i 

PC(NP)| 



0.086 
0.086 
0.09 

0.10-0.32 
0.12 

0.11-0.22 

0.10 

0.10 



1.4 
2.0 
1.5 

2.4- 3.8 
1.0-1.9 

2.5- 2.6 
1.6 

1.2-2.0 



0.69-0.83 

0.55, 0.7 

0.57-0.70 

0.3-0.8 

0.67-0.86 

0.55-0.81 

0.58-0.83 

0.60-0.80 



Nf = 3 
SKS(I)P 



0.14 



0.50-0.94 



this paper, the following abbreviations are used 
to distinguish different lattice actions in simula- 
tions. The first letter stands for the gauge ac- 
tion: P — Plaquette, R — RG improved(|], S — 
Symanzik improved while the second or more 
letters represent the quark action: W = Wil- 
son, KS = Kogut-Susskind, KS(I) = Improved 
KS@gl C(X) = Clover II with Csw= X, where 
TP = Tadpole estimated value and NP = Non- 
Perturbative value. For example, PW stands for 
the plaquette gauge action and the Wilson quark 
action while RC(TP) for the RG improved gauge 
action and the clover quark action with the tad- 
pole estimated value for Csw- 

Simulations have been carried out at sev- 
eral lattice spacings for the PKSjl^ and the 
RC(TP)|l^ actions. Using data from these sim- 
ulations one can make extrapolations to the con- 
tinuum limit. 

3. Several remarks on unquenched simula- 
tions 

In this section several remarks are given on 
problems associated with unqucnchcd simula- 
tions. 

3.1. Auto-correlations 

Since unquenched simulations are very time- 
consuming, it is better to perform hadron mea- 
surements as often as possible. In such cases, one 
has to multiply statistical errors by ^/2Ti^, where 
Tint is the integrated auto-corrclation time of the 
simulation, or equivalently, one should increase 
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Figure 1. Integrated auto-correlation time for the 



number of iterations, T-^^j"" in the HMC algo- 
rithm, as a function of PS meson mass, for the 
PW action 0,|ri|,|9| and the RC(TP) actionpl. 



the bin-size of the Jack-knife analysis until the 
error reaches a plateau. 

The integrated auto-correlation time for an ob- 
servable O is defined as 



1 



t 



(1) 



where C'^{t) is the auto-correlation function for 
O given by 



C°(t) = {0{s)0{s + t))s - {0{s))s - {0{s + t))s{2) 



with (F{s)), 



Fis). 



The auto-correlation time depends on the ob- 
servable and, of course, the simulation algorithm. 
In unquenched QCD, tt meson propagator at a 
given time, or the number of iterations for inver- 
sions of DW where D is a lattice Dirac opera- 
tor, seems to have a long auto-correlation time. 
In Fig. 1^, Tint for the number of iterations Ninv 
in the HMC algorithm is plotted as a function 
of the PS meson mass, which corresponds to 
mps/my > 0.55, in the case of Wilson-type (Wil- 
son/Clover) quark actions |l^ jl^,|l^,|l^ . We es- 
timate Tj^j"" ~ 5 ~ 40, and similarly obtain 
T^nt — 5 ~ 40. One should keep in mind that 
an accurate estimate of error for auto-correlation 
time requires 0(100 — 1000) x Tint data, so er- 
rors in these figures may well be under-estimated. 
These Tint are less or comparable to r/^j ~ 40 for 
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the KS quark action||T^. It is noted that r^nt^'"" 
for the PW are larger than those for the RC ac- 
tion. This may be caused by the difference of 
gauge actions, lattice spacings or simulation pa- 
rameters in the HMC algorithm. 

From a practical point of view, Tint — 10 is 
short enough, even Tint — 40 is not so bad, since 
recent unquenched simulations can accumulate 
8000 trajectories or more, which correspond to 
roughly 100-400 independent configurations. 

It is expected that r^^^"" increases as the quark 
mass (PS meson mass) decreases, according to the 
scaling formula 

r!^r -A-impsa)-^. (3) 

From data at a ~ 0.09 fm for the PW action in 
the figure, one obtains Z = 1.20(5)[|9), which 
predicts a moderate number, t^I""-' ~ 200, for 
mpg = 140 MeV. Since this number is expected to 
be smaller for the RC action, the auto-correlation 
may not be a major obstruction toward lighter 
quark masses in unquenched simulations. 

It has been observed that the topological 
charge Q has a long auto-correlation time of the 
order of a few hundred trajectories or more in 
unquenched HMC simulations with the KS quark 
action. This makes a measurement of the topo- 
logical susceptibility very difficult. On the other 
hand, recent unquenched simulations give r^'^j ~ 
30 - 40 for the RC(TP) action, at a = O.ll fm 
and mps/wv = 0.55 — 0.81 |l5). This value of 
T^^ is small enough to make the calculation of 
the topological susceptibility possible. This dif- 
ference may be explained by the fact that the KS 
fermion determinant suppresses the change of the 
topological charge at small quark masses due to 
chiral symmetry. Further investigations, however, 
will be needed to confirm this speculation, for ex- 
ample, by employing the domain- wall quark ac- 
tion in unquenched simulations. 

In conclusion the auto-correlation is not a ma- 
jor problem for present unquenched simulations, 
except for the case of the topological charge with 
the KS quark action p]]|. 

3.2. Finite size effects 

While finite size effects on hadron masses are 
known to be much stronger in unquenched QCD 
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Figure 2. Finite size errors on hadron masses 
(mH(i) — mH)/TOH for H =PS,V,N as a func- 
tion of spatial extension L (fm), in the case of 
KS (left)||,|| and clover(right)|l|,0 quark ac- 
tions. 

than in quenched QCD for the KS quark ac- 
tion |22[|2^], recent unquenched simulations |p^ , p^ 
suggest that finite size effects for Wilson-type 
quark actions are milder. 

To compare results from the two types of ac- 
tions directly, hadron masses at a finite extension 
L are fitted by the form 

mii{L)=mii + —, (4) 

which is found to be good for the range of L in the 
current simulations [p2[ , and therefore is used as a 
working hypothesis for the comparison. Relative 
errors, (rn}i{L) — rnH)/™H, for H= PS(Pseudo 
Scalar), V(Vector) and N(Nucleon), are plotted 
as a function of L in Fig. |2| at a ~ 0.09 fm and 
mps/mv — 0.68 for the KS quark action p2| , p3[ 
and at a ~ 0.12 fm and mpg/mv — 0.72 for 
the clover quark action[|lj,0. At similar lattice 
spacings and mps/my, the size effect seems to be 
larger for the KS quark. 

Using data currently available, the lattice size 
L (fm) which gives a 3% relative finite size er- 
ror, {mii(L) — mii)/mii ~ 0.03, is summarized in 
Table |. 

From Fig.|| and table |[ one observes that the 
size effect is larger for the KS quark action than 
for the Wilson-type action. Moreover patterns 
of the size effect are different: size effects of PS 
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Table 2 

Lattice size L(fm) which give 3% relative errors. 



mps/my a(fm) 


PS 


V 


N 


KS 








0.5 0.2 


1.2 


1.9 


2.2 


0.09 


2.6 


2.2 


2.6 


0.68 0.09 


2.2 


2.0 


2.3 


Wilson-type 








0.72 0.12 


1.4 


1.5 


1.8 


0.75 0.12 


1.0 


1.4 




0.8 0.12 


1.1 


1.3 


1.7 



and N are comparable and larger than that of V 
for the KS quark action while the size effect of N 
dominates for the Wilson-type quark action. This 
difference is a part of scaling violation since finite 
size effect is universal in the continuum limit, and 
it may be caused by the chiral properties of the 
two actions: the KS quark retains chiral symme- 
try while the Wilson-type quark breaks it explic- 
itly. Another possibility is that the "effective" 
lattice size for the KS quark may be smaller than 
the lattice size for the Wilson quark, due to the 
spread of the KS quark over a hyper-cube in the 
spinor-flavor interpretation. Obvious disadvan- 
tages of Wilson-type quarks such as the absence 
of chiral symmetry, 0{a) scaling violation and 
larger costs of numerical simulations are partly 
compensated by this milder finite size effect. 

3.3. Common scale 

One has to fix a scale to compare results in 
different simulations. In quenched simulations, 
it may be easier and better to use a scale de- 
termined only by the gauge action such as the 
string tension. In unquenched simulations, how- 
ever, all quantities depend on quark actions and 
quark masses as well as the gauge action and the 
gauge coupling constant. 

For the common scale in unquenched simula- 
tions, one had better to use the Sommer scale 
rop4], which is given by 



,dV{r), 



1.65, 



(5) 



where V{r) is the static quark potential. 

This quantity, estimated as ro ~ 0.49 fm, is 
insensitive to either the short-distance physics of 
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Figure 3. (Left) The quark mass dependence 
of the inverse of the Sommer scale for the PKS 
action|p5). (Rightj_The same for PC(TP)0, 
PW|§7"^RC(NP)0 and RC(TP)jl| actions. 
Lines are fits as explained in the text. 



QCD such as the Coulomb term or the long dis- 
tance one represented by string breaking. In sim- 
ulations To is accurately determined by numerical 
derivative of the potential[Q or by the relation 
Tq = (c-l- 1.65)/cr, where the string tension a and 
the Coulomb coefficient c are extracted from the 
fit of V{r). Two methods give consistent values 
of ro with reasonably small errors [T^. 

The quark mass dependence of the Sommer 
scale is shown in Fig.|^, where the vertical axis 
is the inverse of the Sommer scale, while the hor- 
izontal one is the bare quark mass for the KS 
quark action [p5[ or the PS meson mass squared 
for Wilson- type quark actions ||l|,|l|,|l^,|6). One 
can easily obtain the value of vq in the chiral limit, 
fitting the data for small enough quark masses by 
the linear form 



_ + A • M, 

Ro ro 



(6) 



iriaa or m%aa? . 



where M 

Later the Sommer scale at M 7^ becomes nec- 
essary. For distinction I use Rq in such cases in- 
stead of ro, which stands for exclusively the value 
in the chiral limit. 
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3.4. Changing quark mass 

In unquenched simulations, it is non-trivial to 
change the quark mass while keeping the lat- 
tice spacing constant. In a standard method one 
varies the bare quark mass m^a at fixed /3||lO|- 
^0 , determining the lattice spacing in the chi- 
ral limit or at the physical point of mps/my = 
0.135/0.77. In recent simulations (l|Jl| one em- 
ploys a new method in which both niqa and (3 
are varied to keep i?o constant, claiming that the 
lattice spacing becomes smaller at smaller quark 
masses otherwise. In this subsection I will discuss 
if this claim is true. 

First of all it is noted that the lattice spacing 
at different quark masses can not be compared in 
principle since no experimental input is available 
away from the physical quark mass. 

Let me discuss a relation between a and bare 
parameters in QCD using the running coupling 
constant g^{L) instead of Rq to see a difference 
of two methods, since the perturbative expansion 
in an 0(a) improved theory is available for g'^{L) 



as|27 



gl+g^[f3o\og{L/a) + f{mRL) 

+ k mRa\+0{a\gl), (7) 



where g^ is the bare coupling constant, L is the 
lattice size, uir is the renormalized quark mass. 

The mass-independent renormalization scheme 
is defined through the renormalized coupling con- 
stant gj^, defined by 



(8) 



where /i is the renormalization scale. This scheme 
can be easily related to the MS scheme. Clearly 
the standard method of unquenched simulations 
corresponds to this scheme, since a depends only 
on gQ for a fixed g^. The running coupling be- 
comes 



fiL) 



+ k mRa] + 0{a^,g''R), (9) 



which has m^a scaling violation even in the 0(a) 
improved theory, in addition to the physical mass 
dependence /(mijL). Therefore the analogy sug- 
gests that i?o has scaling violation of rriRa as well 



as the physical quark mass dependence in this 
scheme. 

The mass-dependent scheme is specified by 

ffl = g^[l + gU- Po^osiH + firriRL) 

+kmga}]+0{g^„), (10) 

where is the bare quark mass. The lattice 
spacing a depends on both gg and m^a for a fixed 
in such a way that the running coupling be- 
comes mass- independent: 

f{L) =gl + g%pML^i) + 0{a',g%). (11) 

Up to O(a^) ambiguity one can define the rela- 
tion between g^ and m^a so that g^{L) instead 
of g|j becomes constant. In other words the new 
method of unquenched simulations corresponds 
to this renormalization scheme, which keeps 0(a) 
improvement in g^ while has a complicated rela- 
tion to MS scheme. 

This consideration concludes that two meth- 
ods of changing quark mass in unquenched sim- 
ulations correspond to different renormalization 
schemes and that the lattice spacing is kept con- 
stant in both methods, contrary to the claim of 

has proposed another renormalization 
called an improved mass-independent 



rcf.lTe 




Ref. 
scheme 

scheme, whose renormalized coupling constant is 
defined by 



9r = .90 [1 - g'oPo log(Ma)] + Oig^' 



(12) 



where gg — (1 + bgmqa)gQ is called an improved 
bare coupling constant with bg = kg^ -\- 0((7o). 
In this scheme a depends on both g^ and m^a 
through g^ for a fixed while the relation of 
this scheme to the MS scheme remains simple 
since the mass-dependence in vanishes as Oia). 
The running coupling constant in this scheme has 
a mass-dependence and is 0(a) improved: 



f{L) 



~9l- 



-~gy^(3o\og{Lti) + f{7nRL)] + ---.{Vi) 



Let me close this section by suggesting two 
options for unquenched simulations in the fu- 
ture. One can keep using the standard mass- 
independent scheme with the special care about 
m^a scaling violation, which may not be so large 
for light and strange quarks. One can instead try 
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the improved mass-independent scheme with ei- 
ther the perturbative estimate for bg, where the 1- 
loop value is available [2^ or the non-perturbative 
estimate for b^,, which seems possible by the 
method of ref . . 

4. Hadron spectra: Wilson vs. KS quarks 

In order to disentangle the dynamical quark ef- 
fect from the scaling violation, one has to com- 
pare unquenched results with quenched ones in 
the continuum limit. Although precise quenched 
hadron spectra have already been obtained with 
the PW action j^], it is better to confirm the re- 
sults with different actions. Such a comparison 
becomes more important in unquenched QCD, 
where it is much harder to control the continuum 
extrapolation. Therefore, in this section, hadron 
spectra in the continuum limit are systematically 
compared between Wilson and KS quark actions. 

4.1. Edinburgh plot in quenched QCD 

In the continuum limit of quenched QCD the 
nucleon mass is found to be a little smaller than 
the experimental value for the PW action||l[ while 
it is consistent with the experimental value for 
the KS action^^. To see if this small difference 
is caused by the ambiguity associated with the 
chiral extrapolation, the Edinburgh plot(mN/™v 
vs. mps/my) in the continuum limit of quenched 
QCD is shown in the upper part of Fig. |^, for 
PKS@, PW§ and RC(TP)|| actions. 

Surprisingly results between Wilson and KS 
quark actions in the continuum limit of quenched 
QCD do not agree at all in the whole range of 
quark masses. It should be noted, however, that 
the Edinburgh plot from the PW action is con- 
sistent with the one from the RC(TP) action, en- 
suring an agreement among Wilson-type actions 
in the continuum limit. 

Neither the finite size effect nor the chiral ex- 
trapolation can be a main source for this discrep- 
ancy, since the volume is large enough {La ~ 3 fm 
for PW, 2.6 ~ 3 fm for PKS) and the discrepancy 
is significant in the whole range of quark masses. 

To investigate the quality of the continuum ex- 
trapolation, ton/jtiv is shown in the lower part 
of Fig. H as a function of the lattice spacing mya 
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Figure 4. (Upper) The Edinburgh plot(mN/TOv 
vs mps/mv) in quenched QCD at a = 
for PKS|0|, PW|ll and RC(TP)|l|] actions. 
(Lower) m^/rriY as a function of mya at 
mps/mv = 0.5 for same actions. 

at mps/my — 0.5, denoted by the vertical line in 
the Edinburgh plot, where both finite size effect 
and ambiguity of the chiral extrapolation should 
be small. 

Although the quenched PW result has larger 
scaling violation than the quenched RC(TP) re- 
sult, they agree in the continuum limit, conclud- 
ing that the continuum extrapolation is controlled 
for Wilson-type quark actions. The scaling vio- 
lation of the quenched PKS result is well fitted 
quadratically in mya, which is a leading behav- 
ior theoretically expected for the KS quark ac- 
tion||3l[]. It is therefore clear that the discrep- 
ancy between quenched results is not caused by 
the continuum extrapolations. I have to conclude 
that the discrepancy remains unexplained. 

4.2. Edinburgh plot in unquenched QCD 

The analysis in the previous subsection is ap- 
plied also to unquenched results for PKS||l^ and 
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RC(TP)|T5| actions, where the finite size effect is 
expected to be significant, in particular for Ughter 
quark masses. The upper part of Fig. || shows 
the Edinburgh plot in the continuum limit, which 
gradually flattens for the RC(TP)|l|] while con- 
tinuously decreases for the PKS||l^ toward lighter 
quark masses, together with results directly ob- 
tained by simulations at non-zero a without ex- 
trapolations. The behavior for the RC(TP) indi- 
cates the existence of the finite size effect as well 
as the difficulty of the chiral extrapolation from 
the restricted range of quark masses such that 
mps/mv > 0.55. On the other hand, the behav- 
ior for the PKS, which is irregular at non-zero, in 
particular larger a, makes both chiral and contin- 
uum extrapolations difficult. Consequently, the 
quality of data at lighter quark masses for both 
actions is insufficient for the reliable continuum 
extrapolation in the Edinburgh plot. Therefore I 
concentrate on results at heavy masses, where the 
disagreement between the two actions is already 
manifest. 

The lower part of Fig. |^ shows mj^/my as a 
function of the lattice spacing mya at wps/wv = 
0.5, which is obtained by the interpolation for 
the PKS action or by a little extrapolation for 
the RC(TP) action. The scaling violation is 
not so significant for both RC(TP) and PKS ac- 
tions. In particular the result for the latter is 
well fitted quadratically in mya. Two results 
largely deviate even at non-zero a, and the dif- 
ference stays almost constant toward the contin- 
uum limit. Therefore the deviation is not brought 
into by the continuum extrapolation, though the 
linear continuum extrapolation for PC(TP) data 
might be affected a little by higher order terms 
such as a log a or a^. The lattice size. La ~ 2.5 
fm for the RC(TP) or 2.4 fm for the PKS, may not 
be large enough at mps/my — 0.5 in unquenched 
QCD. As mentioned in Sect. p.l] , there may be 3% 
error in ton for the KS quark action at La ~ 2.6 
fm and a ~ 0.1 fm. Although no data are avail- 
able for the Wilson quark action at this mps/my, 
the error is expected to be smaller than the KS 
quark action. With all systematics together, it 
may be too early to conclude that results between 
Wilson and KS quark actions disagree in the con- 
tinmmr limit at mps/my = 0.5. Further inves- 
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Figure 5. (Upper) The Edinburgh plot in un- 
quenched QCD for the PKS action(diamonds) [p^ 
and the RC(TP) action(circles)[^. Open sym- 
bols represent data at a = 0.32, 0.28, 0.23, 0.15, 
0.10 fm for the former and 0.22, 0.16, 0.11 fm 
for the latter, from above to below, while solid 
symbols correspond to results in the continuum 
limit. (Lower) m^/my as a function of my a at 
wps/wy = 0.5 for same actions. 

tigations, in particular on the finite size effect or 
at larger mass ratios such that nips/my = 0.6 
to avoid the chiral extrapolation, are urgently 
needed to draw definite conclusions. 

4.3. Hadron spectra normalized by Rq 

In the previous subsection, the discrepancy in 
the Edinburgh plot between Wilson and KS quark 
actions became manifest in the continuum limit, 
in particular, for the quenched case. However, 
since the Edinburgh plot represents the relation 
among mass ratios, it is difficult to find a source 
responsible for the discrepancy. Alternatively, 
in this subsection, hadron masses normalized by 
the Sommer scale are considered as a function of 
(mpsi?o)^, which plays the role of quark mass. 
With this normalization one can keep the sim- 
ple relation between toh and long as 
the quark mass dependence of Rq is sub-leading. 
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(mpsRo) 



Figure 6. myRo as a function of (mpsi?o)^ in un- 
quenched QCD for the PKS action(diamonds) 
and the RC(TP) action(circles)|l|]. Open sym- 
bols represent data at non-zero a, while solid 
ones in the continuum limit. Solid line represents 
quenched data in the continuum limit for the PW 
action [0 for comparison. 

while still avoiding the chiral extrapolation. 

In Fig. 1^ myRo from unquenched simulations 
is plotted as a function of (rnpsi?o}f at non-zero 
a, for both RC(TP)[|5| and PKS@ actions, to- 
gether with lines in the continuum limit from 
quenched PW|| and unquenched RC(TP). Un- 
fortunately lines in the continuum limit are not 
available for the PKS action. 

Data at non-zero a approach the continuum 
limit from below(a = 0.22, 0.16, 0.11 fm) for 
RC(TP) and from above(a = 0.21, 0.15, 0.11 fm) 
for PKS. The behavior looks reasonable except 
for PKS at a ~ 0.21 fm. Roughly speaking slopes 
are similar among unquenched results and steeper 
than the quenched one. 

As before the scaling behavior is examined at 
(mps-Ro)^ = l-22, where PKS data are available. 
In the left part of Fig. ^my-Ro at (mps-Ro)^=l-22 
is shown as a function of u/Rq. The quenched 
PKS datag and SKS(I) data||,|| agree with 
each other in the quadratic continuum extrapo- 
lation, and an extremely good scaling behavior 
is observed for the latter. The continuum values 
also agree with that for the PW action. There- 
fore no discrepancy is found for quenched vector 
meson masses between KS and Wilson quarks. 

Similarly the scaling behavior of ton^o 
is shown in the right part of Fig. |^ at 
(mpsi?o)^ = l-22. The linear continuum extrap- 
olation of quenched PW data||l[] gives a smaller 
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Figure 7. mvi?o(left) and mNi?o (right) as 
a function of a/Rg at {mpsRo}^ = 1.22, 
in quenched QCD(open symbols) with PW|0], 
PKS|^ and SKS(I)|2|,|| actions and in un- 
quenched QCD(solid symbols) with RC(TP)||l|] 
and PKSOl actions. 



value than the quadratic one of quenched PKS 
which well agrees with the one of SKS(I) 



data|30 



data|32 



The individual analysis for my and ttin sug- 
gests that a discrepancy in ton is a main source 
for the disagreement in the Edinburgh plot be- 
tween PW and PKS actions in quenched QCD. 

For unquenched QCD, continuum extrapola- 
tions of my and ton are reasonably good for 
the RC(TP) actionjlSj and agree with those for 
quenched QCD with PWljl, while very large scal- 
ing violations for my and ton , which do not show 
expected behaviors, make continuum extrap- 
olations of individual masses impossible for the 
PKS action[|l3|. However these large scaling vi- 
olations in my and ton cancel each other in the 
ratio ton/tov as seen in the lower part of Fig. ^. 

4.4. Suggestions on further investigations 

From available data so far we can not iden- 
tify a definite reason for the discrepancy in the 
Edinburgh plot between Wilson and KS quark 
actions, in particular for quenched QCD, where 
the disagreement is found only in ton. To com- 
pletely resolve the problem, we should accumu- 
late as many data as available for all varieties of 



9 



actions and simulations. 

We should analyze the Edinburgh plot and the 
plot of itihRo vs (mpsi?o)^ , as well as their scaling 
behavior at mps/mv=0.5, 0.7 for the former and 
at {mpsRo)^ — 1-22, 3.0 for the latter as typical 
values already used in several literatures. 

Even quenched results from the third type of 
quark actions such as domain- wall or overlap will 
become useful to resolve the problem. 

5. Dynamical quark effects on hadron 
spectra 

In this section, dynamical quark effects on 
hadron spectra are investigated, comparing data 
from unquenched simulations with those from 
quenched ones. Due to the discrepancy between 
Wilson and KS quark actions in the previous sec- 
tion, results only for Wilson-type quark actions 
are considered here, since more data are available 
so far. 

5.1. Quark mass dependence in partially 
quenched analysis 

In unquenched simulations, valence quark 
masses can be different from the mass of dynam- 
ical quarks (sea quark mass). In such a partially 
quenched analysis, the PS meson mass squared 
are employed for Wilson-type quark actions to 
specify how large the quark mass is, without chi- 
ral extrapolation necessary to define the quark 
mass through the critical hopping parameter. 

The following equality is then assumed for the 
PS meson mass. 



'PS.val 



(14) 



where s specifies a sea quark mass for degenerate 
Nf = 2 dynamical quarks and vi and V2 denote 
valence quark masses in the meson. This assump- 
tion is correct at the lowest order of chiral per- 
turbation theory and is well-satisfied numerically 
in data of ref. [0 . 

In Fig.^ rnv(s; ''^i: ^^2) is plotted in lattice units 
as a function of mpg ^^^j . At a fixed sea quark 
mass{Ksea — 0.1355, for example, which corre- 
sponds to mps/my ~ 0.6), data lie on one curve 
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Figure 8. mv(s;wi,t'2) in lattice units as a func- 
tion of mpg_^j^i at /3 = 5.2 for the PC(NP) action 
on a 16^ X 48 lattice |l7|. 



with some curvature, suggesting that higher or- 
der terms in m|g ^^^j are necessary while the effect 
of (TOpg(s; wi, wi) — mpg(s; W2, ^12))^ term is negli- 
gible. Therefore, as shown in the figure, data can 
be reasonably fitted by the form 



mv(s;wi,W2) = A + Bs-m 

+ {By + a 



2 

PS,! 



'PS,! 



+ ■ (™ps,val)^ 



™PS,val 

(15) 



with 



■-PS.s 



— mpg(s; s, s), where all possible 
terms of O(mpg) are included. The fact that the 
intercept and the curvature of each curve depend 
on the sea quark mass might be interpreted as a 
dynamical quark effect. 

To see if this is really a dynamical quark effect, 
the same data normalized by Rq are shown in the 
left part of Fig. ||, together with quenched data at 
P = 6.0 and 6.2 for the same action|3^]. Note that 
the scaling violation is negligible in the quenched 
case. 

A curvature can not be clearly seen anymore, 
and a linear decrease seems steeper than the 
quenched one. Partly due to larger errors caused 
by that of Rq, however, the fit by the form 



Rq ■ mv(s;wi,W2) 



A + Bs-Rl- mis 
Bv ■ Rq ■ mpg^^^i 



(16) 



can not resolve a sea quark mass dependence 
in the data, giving Eg ~ 0. The fit with 
Bg ~ 0, plotted in the figure, gives a reason- 
able value of x^, showing a steeper slope _B„, 
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Figure 9. (Left) Romy{s;vi,V2) as a function 
of (TOps,vai-Ro)^ for the PC(NP) action from 
the JLQCD collaboration|p7|] , together with the 
quenched data for the PC(NP) actionQ. Sym- 
bols are the same in Fig. ^. (Right) The same plot 
for the PC(NP) action from the UKQCD collab- 
oration together with the fit line of the left 
figure. 

which may be regarded as the dynamical quark 
effect. Further precise calculations will be re- 
quired, however, to resolve the sea quark mass 
dependence of the slope Bs in this range of sea 
quark masses(mps/TOv — 0.6 — 0.8), so that the 
dynamical quark effect is undoubtedly identified. 

As mentioned before, the UKQCD collabora- 
tion has performed an unquenched simulation 
varying both /3 and m^a of PC(NP) action in or- 
der to keep Rq constant [|6|. As shown in the 
right part of Fig. ||, the result is consistent with 
the one from the fixed f3 simulation of the JLQCD 
collaborationJlTt , once both data are normalized 
by Ro- 

5.2. Chiral extrapolation 

The dynamical quark effect is expected to 
be enhanced toward smaller sea quark masses 
reached by the chiral extrapolation, which is 
briefly described in this subsection. 

The light quark mass is determined by the 
experimental value of mjr/mp. In the present 
method, the ratio determines the corresponding 
PS meson mass, wps.sO by solving the equation 

_ rnps.sfl 

rup ~~ A+ {Bs + B^)m1,ga^ H 



where the denominator is obtained by setting 
™PS vai = ™PS,s dll)' ^^fi lattice spac- 

ing a is explicitly written. Hereafter mpg ^a^ is 
fixed to the solution of this equation, mps,sa = 

There are several ways to determine the strange 
quark mass. The if-input uses the experimen- 
tal value of mif/TOTT to determine other strange 
mesons, while the 0-input uses the experimental 
value of m^/mp 

To compare results at different a, hadron 
masses are normalized by tq, the value of i?o in 
the chiral limit. 

5.3. Meson spectra at dynamical light 
quark masses 

The vector meson mass normalized by ro, 
myro in unquenched QCD is shown in the left 
part of Fig.|l^, as a function of rnpg^^^rQ, at 
a = 0.11 - 0.22 fm for the RC(TP) action|l|l 
and at a = 0.1 for the PC(NP) action||l|l, to- 
gether with the quenched result in the contin- 
uum limit 1^. The lines in the figure are de- 
termined by eq. (^ with mpg ^ = fixed, 
and the points on the line marked by (7r,p), 
{K,K*) and {ris,4>) are given by mps/my = 
0.1350/0.7684, 0.4977/0.8961 and 0.695/1.0194, 
respectively, where m^^ =0.695 GeV is used for an 
unphysical rj^. The error of the line, not shown in 
the figure for the RC(TP) action, is roughly equal 
to the error of the points, where the error of tq 
is not included. The error band for the PC(NP) 
action includes this source of error. 

The slope of unquenched data increases as a 
decreases, so that the slope at a = 0.11 fm of 
RC(TP) action is clearly steeper than the slope 
of the quenched result in the continuum limit, and 
closer to experiment. This is the manifestation of 
dynamical quark effect on hadron spectra we are 
looking for. The chiral extrapolation in the sea 
quark mass, as well as control of scaling violation, 
is necessary to see the effect clearly. The line for 
the PC(NP) at a = 0.1 fm, which lies a httle 
above, seems consistent with the RC(TP) if errors 
from the remaining scaling violation and rg are 
taken into account. 

The result of my can be converted into the 
hyper-fine splitting of mesons, and is shown in 
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Figure 10. (Left) TOv?'o(left) and (Right) 
(my — mpg ^^ifT^ as a function of mpg vai'^o with 
"^PS s ~ ™J fixed, in unquenched QCD for 
RC(TP)(l| and PC(NP)|3 actions, together 
with the quenched result of the PW action Q. 



the right part of Fig. |T^, where (my 
is plotted as a function of ^j^\T^ with mpg ^ — 
fixed. The slope of unquenched data de- 
creases with a, and is very different from the 
quenched one in the continuum limit. The dy- 
namical quark effect is most noticeable for the 
hyper- fine splitting normalized by tq. 

Let me remark that a steeper slope of the 
hyper-fine splitting in quenched QCD means a 
smaller hyper-fine splitting of strange mesons 
than the experimental value, if the splitting be- 
tween uip and m^n- is fixed to the experimental 
value. 

5.4. Continuum limit of strange meson 
masses 

In Fig. |ll|, the scaling behavior of and 
K* meson masses from if-input, converted into 
GeV unit using p meson mass, are plotted in 
Nf = 2 unquenched QCD for RC(TP)jl| and 
PC(NP)(l7| actions, together with the quenched 
result for PW@] and RC(TP)|5| actions. Lines 
represent linear continuum extrapolations for the 
Nf = 2 RC(TP), quenched PW and RC(TP) re- 
sults. 

A large discrepancy between the experimental 
value and the quenched result in the continuum 
limit is much reduced in Nf~2 unquenched QCD. 
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Figure 11. and m/f. (GeV) as a function of 
a (GeV) in Nf ^ 2 QCD for RC(TP)[|l| and 
PC(NP)||l|] actions and in quenched QCD for 
PW|] and RC(TP)(l|] actions. The scale is set 
by rup. 

It will be interesting to see if the inclusion of the 
dynamical strange quark removes the remaining 
small difference. 

5.5. Remark on the continuum extrapola- 
tion 

Here a remark is given on the scaling be- 
havior of perturbatively improved theories such 
as RC(TP), extensively employed in unquenched 
QCD simulations. Suppose that a renormalized 
observable X{a) at non-zero lattice spacing a is 
expanded in the bare coupling constant gQ as 

X{a) = Xo+gllx^+aixl+ylloga)] 

+9o [^2 + vl log a + a{x\+yl log a 
+zl\og^a)\+0{a\gl) (18) 

Note that the g^ log a term is absent because of 
the renormalizability. The tree level 0{a) im- 
provement leads to y\ — Q and z\ ^ while 
the 1-loop Ola) improvement gives x\ — Q and 
j/2 = in addition. In terms of the renormalized 
coupling at scale p defined as 



9l 



9lr 



vl 



o5i?.log(/xa) 



9l 



% log(/xa) 
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one obtains 
X{a) 



X{Q) + g\a [x\+y\\oga] 



g%a [x\ + yl \oga- 



4 log^ a] 



The leading terms of scaling violation are g^^^a 
and g'^a log a for unimproved theories such as the 
Wilson fermion, g|j.a and g'^aloga for tree- level 
or tadpole 0{a) improved theories, or g'^a and 
5^aloga for 1-loop 0(a) improved theories. This 
concludes that a linear continuum extrapolation 
employed for data in the RC(TP) action ||l5| is 
justified if the g\a term dominates the (?|.a log a 
term and O(a^) terms can be neglected. 

Note that terms such as (gn)"a ~ t, ^ 

^ (log a)" 

never appear in the scaling violation, since the 
physical scale /i is independent on a. 

5.6. Dynamical quark effect on baryon 
masses 

Finally the dynamical quark effect on baryon 
masses is briefly discussed. The chiral extrapola- 
tion of baryon masses is made in terms of the PS 
meson mass, by the formula explicitly given as 

mY.{s;vi,V2,V2) = A + Bg^s + 2F^y^ 
+ (F- D)fiy, + Csfil + {E + Eo)fil, 



+ {E - Eo)^J.l, + iG + Go)^J-s^J-v^ 

+ {G — Go)f^sf^V2 + Gy^2 



(19) 



for S-like octet baryons, where ^.s — Wpg g and 
{s;vi,Vi). Similar formula are given 
for the A-like octet baryon mass, mA(s; vi,V2,V2), 
and simpler ones for the decuplet baryon mass 
mDis;vi,v2,V3)^. 

In Fig. ^ the E-like octet baryon masses are 
plotted as a function of the averaged PS meson 
mass squared, ?7ips vai ~ (/^"i + 2/1^,3) /3, for the 
RC(TP) action|l|.' 

To draw lines in the figure, the dynamical 
quark mass is kept fixed to the light quark 
mass, fis = rn^. The solid line corresponds to 
roTOE(s; w, D, d) as a function of r§yUt,, while the 
solid symbols represent N at the point oi v = s. 
The error of each line, which does not include 
the error of ro , is roughly equal to the size of the 
symbol on the line. 




PS,val 



Figure 12. S-like octet baryon masses, normal- 
ized by ro, as a function of tJ mpg ^^^j for the un- 
quenched RC(TP) action[^ and the quenched 
PW action §. 



The scaling violation is small at a < 0.16 fm 
and the slope of lines in unquenched QCD be- 
comes flatter than that in quenched QCD, as a 
decreases. This may be a dynamical quark ef- 
fect. One must, however, be careful to conclude 
this, since the physical lattice size in the simula- 
tions are different, ~ 2.5 fm in unquenched sim- 
ulations and 3 fm in quenched ones, so that fi- 
nite size effects could make such a difference. It 
is therefore interesting to obtain the correspond- 
ing plot in quenched QCD with the lattice size 
of 2.5 fm in the continuum limit, to identify the 
dynamical quark effect on the slope. Moreover, 
even though some finite size effect is expected to 
exist on baryon spectra, the small scaling viola- 
tion enables us to predict the chiral behavior of 
the baryon masses in the continuum limit with a 
fixed lattice size of 2.5 fm. 

Similar results are obtained for the A-like octet 
baryon mass and the decuplet baryon mass. In 
particular almost no scaling violation is observed 
for the decuplet baryon mass at a < 1.6 fm. 

6. Dynamical quark effect on other quan- 
tities 

6.1. String breaking 

In the past unquenched simulations, the ex- 
pected behavior for string breaking that the static 
quark potential V{R) levels off at large enough i?, 
has not been observed] 35 -3^: V{R) increases lin- 
early in R beyond the point where it exceeds the 
energy of the two mesons, due to poor overlap of 
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Figure 13. Lowest 3 energy levels, corresponding 
two S-wave mesons, S+P wave mesons and the 
string, as a function of the distance R. 

Wilson loops with the meson-meson state. 

Several variational methods have been intro- 
duced to find better operators which have larger 
overlap with the meson-meson state, and success- 
fully applied to several models such as an adjoint 
Wilson loop in pure gauge theories or gauge-Higgs 
models, showing the expected behavior for "string 
breaking" ||. 

At lattice'99 the MILC collaboration has re- 
ported a first evidence of string breaking in un- 
quenched QCD simulations! 38 1, and the follow-up 
of the result is presented at this conference p9[. 
They employ the Nf = 2 KS fermion for dy- 
namical quarks on a (3.3fm)'^ x (3.9fm) lattice at 
a ~ 0.16 fm and mps/wv = 0.36. The energy ob- 
tained from meson-meson as well as meson-string 
correlators starts to level off at i? ~ 1 fm, where 
the energy is roughly equal to twice the mass 
of the heavy- light meson|38|. The lowest three 
eigen-values, obtained by the diagonalization of 
correlators including the string-string one, are 
shown in Fig. |l^|^, which supports the result 
in ref. iQ. The mixing between the string state 
and the meson-meson state, however, is found to 
be too weak to show level crossing 

Let me give one comment on the definition of 
string breaking. String breaking should be a phe- 
nomenon that the energy for one i?-independent 
operator linearly increases at short distance and 



levels off at long distance. According to this cri- 
terion the behavior in ref. is really string 
breaking. If one uses different operators at dif- 
ferent R instead, the lowest energy can always 
be obtained by a diagonalization, even in the ab- 
sence of the transition from the string state to the 
meson-meson state. 

Apart from the string breaking, decays of 
hadrons, such as p ^ tttt, give some indications of 
the dynamical quark effect, since such decays are 
possible but largely suppressed in the quenched 
approximation. At this conference, the MILC col- 
laboration presents a result which indicates the 
decay of ao(0++) -> 7r(0-+)77(0-+)@. Contrary 
to the case of string breaking, the operator 
employed in the simulation seems to well overlap 
with both oo and nr] states. 

6.2. Topological susceptibility 

From the flavor-singlet 

axial Ward-Takahashi(WT) identity in the con- 
tinuum. 



{id>'A^{x) + 2m,P{x))Q{y)) 
2Nf{Q{x)Q{y)) = 0, 



one obtains 



2mq / d^x {P[x)Q) 



2Nf{Q^ 



(20) 



(21) 



where Q = J d^x Q{x) is the topological charge 
with the topological charge density Q{x) = 

^\f^Ax)F^^''{x). Therefore 



V 



NfV 



d'^x{P(x)Q) ^ (22) 



as niq 0. The topological susceptibility x must 
vanish in the chiral limit. An investigation of 
this property on the lattice becomes a good check 
for lattice chiral symmetry in unquenched simula- 
tions. The decrease of x toward mg = is driven 
purely by the fermion determinant. 

In the past simulations, the expected decrease 
of X has not been observed for both KS[^ and 
Wilson- type fisll quark actions. The topological 
susceptibility normalized by Rq does not show the 
quark mass dependence at all and values are con- 
sistent with the quenched onellS]. 
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Figure 14. xRq as a function of rupg gi?g in Nf ~ 
2 unquenched QCD for the PC(NP) action(open 
circles) [H at a ~ 0.1 fm and for the RC(TP) 
action(sohd symbols) at a ~ 0.11 fm(circles)||4^ 
and 0.16 fni(squares) Jisf . 

This year the UKQCD collaboration^] and 
CP-PACS collaboration have obtained new re- 
sults, which are shown in Fig. |lj. The decrease 
of X toward the chiral limit is observed at a ~ 0.1 
fm for both PC(NP) and RC(TP) actions, in con- 
trast to the previous result at a = 0.16 fm for the 
RC(TP) action, though it is still far from verify- 
ing the flavor-singlet axial WT identity. Theoreti- 
cal considerations as well as further numerical 
analyses are needed to settle this issue completely. 

6.3. Flavor-singlet PS meson 

The fact that the flavor-singlet pseudo-scalar 
meson is not a Nambu-Goldstone boson(U(l) 
problem) must be demonstrated by lattice QCD 
simulations. Numerically this has been a diffi- 
cult task since one has to calculate the quark 2- 
loop diagram special for the flavor-singlet meson. 
By the volume-sources method without gauge- 
fixing^^ or the U(l) random noise method fl^ 
the estimation of the flavor-singlet meson mass 
becomes possible in quenched QCD. 

The calculation of the singlet meson mass ap- 
pears to be much harder in unquenched QCD, 
where the individual meson propagators should 
behave as 



M0)vit))2-\oop ^ Z.e-™'* - Z.e-™"* (23) 
{mvit))i-loov ^ ^.e-™-*. (24) 
Here rj — q^slq is the flavor-singlet PS meson 
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Figure 15. to,, in GeV as a function of a in Nf — 2 
unquenched QCD for RC(TP)|l|, PWH and 
PKS^^. The smearing source is employed for 
the latter two. The horizontal line is the "experi- 
mental" value of the flavor-singlet meson mass in 
Nf = 2 QCD. 

operator, is the corresponding mass and to^ is 
the mass of the degenerate non-singlet PS meson. 
Since the pole of the appears after the leading 
contribution of tt is completely canceled between 
two diagrams, the signal becomes too noisy at 
large t to extract m^. 

Therefore one is forced to extract m,, at small 
t, and one usually fits the ratio of propagators 

asira 



(^(0)^W)2-loop 
('?(0)'?(i)>l-loop 



1 



-niT^)t 



hoping that a large part of contributions from ex- 
cited states may be canceled out in the ratio. The 
CP-PACS collaboration[|l5| indeed found that 
rriri — rriTr extracted from the ratio for tmin = 2 is 
consistent with the one for tmin = 3, 4 within sta- 
tistical errors, where tmin is the minimum value 
of t for the fit . 

Recently the SESAM-T^L collaboration has 
shown I^tI that the plateau for to,, as well as to^ 
starts at reasonably small i, employing the smear- 
ing source. This opens a possibility for a reliable 
calculation of in very near future. 

In the chiral limit one obtains a non-zero value 
of to,, |^ , P^ , ^ , which is shown in Fig. [Tsl as a 
function of a, together with the value of to,, ~ 0.7 
GeV for Nf = 2 QCD derived from the experi- 
mental value(iV/ = 3) by the Witten-Veneziano 
formulajl^. Although results from different ac- 
tions seems consistent within large errors, the re- 
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suit from the RC(TP) action show a systematic 
increase toward the continuum Hmit. One has 
to check if this large scaling violation is caused 
by the contamination from excited states to the 
ground state in the ratio R{t). 

A new calculation with the smearing source on 
configurations of the CP-PACS collaboration will 
make it possible to extract reliably in the con- 
tinuum limit of Nf = 2 QCD, while the inclusion 
of the mixing to ss state as well as the dy- 
namical strange quark will be necessary for the 
final result in iV/ = 3 QCD. 

7. Summary 

We surprisingly find that hadron spectra from 
quenched QCD disagree in the continuum limit 
between KS and Wilson quark actions. This dis- 
crepancy has to be resolved as soon as possible. 
Personally I feel that some problems may exist 
in the identification of particles in KS quark ac- 
tions, in particular of baryons. The restoration 
of the flavor breaking in hadron spectra might be 
a good check. Hadron spectra from domain-wall 
and overlap quark actions |Q may also give an 
important hint for the solution. 

Within Wilson-type quark actions, the Nf — 2 
dynamical quark effect on hadron spectra, in par- 
ticular in meson masses, has been demonstrated. 
Both chiral and continuum extrapolations are im- 
portant to identify the effect. Results for baryons 
suggest that the scaling violation is small but the 
lattice size has to be increased for reliable pre- 
dictions in the continuum limit. Decreasing the 
dynamical quark mass will also be the next chal- 
lenge in unquenched simulations. 

The dynamical quark effect becomes manifest 
in some other quantities such as the string break- 
ing, the topological susceptibility and the flavor- 
singlet meson mass, even at finite lattice spacing. 

New direction of investigations in unquenched 
QCD will be opened by realistic Nf — 2+1 un- 
quenched simulation, which is the last step to 
numerically solve QCD. Promising algorithm for 
Nf = 2 + 1 QCD seems HMC for Nf = 2 hght 
quarks and PHMC||5|] for the Nf = 1 strange 
quark. Some results will be expected at the first 
lattice conference in the next century. 



Acknowledgements 

I would like to thank C. AUton, J.M. Carmona, 
A. Hart, A. Irving, K. Ishikawa, D.B. Leinweber, 
J.W. Negele, H. Panagopoulos, R. Parthasarathy, 
G. Rossi, S. Sharpe, D. Toussaint, U. Wolff and T. 
Yoshie for communicating and providing their re- 
sults, and, in particular, R. Burkhalter, C. DeTar, 
S. Gottlieb and T. Kaneko in addition for mak- 
ing new data I have requested. I also thank R. 
Burkhalter, Y. Iwasaki, T. Kaneko, S. Hashimoto 
and A. Ukawa for valuable comments on the 
manuscript. This work is supported in part by 
the Grants-in-Aid of Ministry of Education(Nos. 
12014202, 12640253). 

REFERENCES 

1. S. Aoki et. al. (CP-PACS collaboration), 
Phys. Rev. Lett 84 (2000) 238. 

2. V. Lubicz, "Quark mass on lattice: Light and 
heavy" , these proceedings. 

3. L. Lellouch, "Light hadron weak matrix ele- 
ments" , these proceedings. 

4. C. Bernard, "Properties of heavy-light and 
heavy- heavy hadrons, and constraints on the 
unitarity triangle", these proceedings. 

5. Y. Iwasaki, Nucl. Phys. B258 (1985) 141; 
Univ. of Tsukuba report UTHEP-118(1983), 
unpublished. 

6. M. Liischer and P. Weisz, Commun. Math. 
Phys. 97 (1985) 59; P. Weisz, Nucl. Phys. 
B212 (1983) 1. 

7. S. Naik, Nucl. Phys. B316 (1989) 238. 

8. K. Orginos, D. Toussaint and R.L. Sugar, 
Phys. Rev. D60 (1999) 054503; Nucl. Phys. 
B(Proc. Suppl.) 83-84 (2000) 878. 

9. B. Shcikholcslami and R. Wohlert, 
Nucl. Phys. B259 (1985) 572. 

10. N. Eicker et. al. (SESAM collaboration), 
Phys. Rev. D59 (1999) 014509. 

11. T. Lippert et. al. (SESAM-T^L collabora- 
tion), Nucl. Phys. B(Proc.Suppl.) 60A (1998) 
311. 

12. R. Mawhinney, Nucl. Phys. B(Proc. Suppl.) 
83-84 (2000) 57 and references therein. 

13. C. Bernard et. al. (MILC collaboration), 
Nucl. Phys. B(Proc.Suppl.) 73 (1999) 198; 



16 



14. 



15. 



16. 



17. 



18. 



19. 



20. 



21, 



60A (1998) 297 and references therein. 
C.R. AUton et. al. (UKQCD Collaboration), 
Phys. Rev. D60 (1999) 034507. 
R. Burkhalter, Nucl. Phys. B(Proc.Suppl.) 73 
(1999) 3; S. Aoki et. al. (CP-PACS collabora- 
tion), ibid. 73 (1999) 192; A. Ali Khan et. al. 
(CP -PACS collaborat ion), ibid. 83-84 (2000) 
176; |hep-lat/0004010| ; these proceedings. 
A.C. Irving (UKQCD Collaboration), these 
proceedings. 

S. Aoki et. al. (JLQCD Collaboration), these 
proceedings. 

C. Bernard et. al. (MILC Collaboration), 
Phys. Rev. D62 (2000) 034503. 
Th. Lippert, S. Giisken and K. Schilling, 
Nucl. Phys. B(Proc.Suppl.) 83-84 (2000) 182. 
Y. Kuramashi et. al, Phys. Lett. B313 (1993) 
425; B. AUes et. al, ibid. B389 (1996) 107. 
The long auto-correlation time of Q with the 
KS action may be reduced by the parallel 
tem pering method. See E.-M. Ilgenfritz et. 
al, |hep-lat/0007"039 . 



22. M. Fukugita et. al, Phys. Rev. Lett. 68 (1992) 
761; Phys. Lett. B294 (1992) 380; Phys. Rev. 
D47 (1993) 4739; S. Aoki et. al., ibid. D50 
(1994) 486. 

C. Bernard et. al. (MILC Collaboration), 
Phys. Rev. D48 (1993) 4419. 
R. Sommer, Nucl. Phys. B411 (1994) 839. 
S. Tamhankar and S. Gottlieb, Nucl. Phys. 
B(Proc.Suppl) 83-84 (2000) 212. 
B. Bolder et. aL(SESAM-TxL CoUabora- 



23. 

24. 
25. 

26. 

tion), |hep-lat/0003012 . 

27. M. Liischer et. aZ. (ALPHA collaboration), 
Nucl. Phys. B478 (1996) 365. 

28. S. Sint and R. Sommer, Nucl. Phys. B465 
(1996) 71. 

29. G. MartineUi et. al, Phys. Lett. B411 (1997) 
141. 

30. C. Bernard et. al. (MILC Collaboration), 
Phys. Rev. Lett. 81 (1998) 3087. 

31. S. Sharpe, Nucl. Phys. B(Proc.Suppl.) 34 
(1994) 403. 

32. C. Bernard et. al. (MILC Collaboration), 
Phys. Rev. D61 (2000) 111502. 

33. C. Bernard et. al. (MILC Collaboration), 
these proceedings. 

34. K.C. Bowler et. al. (UKQCD Collaboration), 



Nucl. Phys. B(Proc.Suppl) 63 (1998) 182; 
Phys. Rev. D62 (2000) 054506. 

35. K. Schilhng, Nucl. Phys. B(Proc.Suppl) 83-84 
(2000) 140 and references therein. 

36. S. Aoki et. al. (CP-PACS collaboration), 
Nucl. Phys. B(Proc.Suppl.) 73 (1999) 216. 

37. B. Bolder et. aL(SES AM-TxL Collabora- 



tion), |licp-lat/0005018 

38. C. Detar, U. Heller and P. Lacock, 
Nucl. Phys. B(Proc.Suppl) 83-84 (2000) 310. 

39. C. Bernard et. al. (MILC Collaboration), 
these proceedings. 

40. P. Penna nen and C. Micha el (UKQCD collab- 
oration), |hep-lat/0001015| 



41. B. AUcs e t. 
lat/0004020| . 

42. A. Hart and M. 



al., |hep-iat/9909045| ; 



Teper (UKQCD collabora- 



tion), |hep-lat/0004180 . 

43. S. Aoki et. al. (CP-PACS collaboration), in 
preparation. 

44. Some subtleties exist in the flavor-singlet lat- 
tice WT-identity: G. Rossi and M. Testa, pri- 
vate communication. 

45. Y. Kuramashi et. al., Phys. Rev. Lett. 72 
(1994) 3448. 

46. G. Kilcup et. al, Nucl. Phys. B(Proc.Suppl) 
47 (1996) 358; L. Venkatamaran et. al, ibid 
53 (1997) 259. 

47. T. Stru ckmann et. aL(SESAM-TxL CoUab o- 
ration), |hep-lat/0006012| ; |hcp-lat/0010005| . 

48. E. Witten, Nucl. Phys. B156 (1979) 269; 
G. Veneziano, ibid B159 (1979) 213. 

49. C. Michael et. aZ.(UKQCD collaboration), 
Nucl. Phys. B(Proc.Suppl.) 83-64 (2000) 162. 
C. McNe ile and C. Michae l (UKQCD collab- 
oration), |hep-lat/0006020[ 

50. P. Vranas, "Domain-wall fermions and appli- 
cations", these proceedings. 

51. T. Takaishi and P. de Forcrand, Nucl. Phys. 
B(Proc.Suppl.) 53 (1997) 968; these proceed- 
ings; R. Frezzoti and K. Jansen, Phys. Lett. 
B402 (1997) 328. 



